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Abstract 

Based on calculations using the local density approximation, we propose 
quantum wire networks with square and plaquette type lattice structures that 
form quantum dot superlattices. These artificial structures are well described 
by the Hubbard model. Numerical analysis reveals a superconducting ground 
state with transition temperatures T c of up to 90 mK for the plaquette, which 
is more than double the value of 40 mK for the square lattice type and is suf- 
ficiently high to allow for the experimental observation of superconductivity. 
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The electronic properties of solids are closely related to their crystal structure which is 
strictly determined by the atomic nature of the individual elements. However, the imposition 
of a superstructure on a given lattice makes the controlled fabrication of structures with 
chosen properties possible. One approach discussed more recently is to place quantum dot 
(QD), also known as artificial atoms [0,0, on the points of a lattice to form an artificial crystal 
called a quantum-dot superlattice (QDSL). Prerequisites for obtaining a band structure 
based on the periodic potential in a quantum-dot lattice are coherent coupling and good 
uniformity in the QD array. Recent experiments have reported the realization of such 
QDSLs. The observation of the quantum Hall effect in a lateral periodic potential based on 
a two-dimensional electron system |3[] is clear evidence of phase coherent transport across 
unit cells. Schedelbeck et al. [|J used cleaved edge overgrowth to fabricate coupled QDs 
that formed in the intersection of three quantum wells. They observed peak-splitting in 
the photoluminescence spectrum demonstrating the coherent coupling of states localized in 
the QDs. The QDs formed because of the larger area at the intersections, which lowers the 
confining potential. In a similar way, QDSLs can be created using two-dimensional quantum 
wire networks which allow the formation of a large number of coupled dots with the high 
regularity of the underlying grid. This type of wire network geometry has already been 
achieved with GaAs H and Si || wires. 

The use of QDSLs should make it possible to observe the electron correlation effects 
predicted in mathematical models of lattice systems, such as the Hubbard model. The 
Hubbard Hamiltonian is defined as H = —J2i,j,atij c ta c jo- + U J2i n ii n i.i, where is the 
hopping integral between different sites i and j, U is the intra-site interaction between 
electrons with opposite spin, and (q) is the electron number (annihilation) operator with 
spin a on site i. Recently, we have used the spin dependent local density approximation 
(LDA) as a basis for proposing a feasible method for achieving a ferromagnetic ground 
state in a Kagome QDSL within a wide electron filling range []7|, which is consistent with a 
mathematical proof || of the ferromagnetism in the Hubbard model on the Kagome lattice. 

In this Letter, we propose a method for forming QDSLs in a quantum wire network of 
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square and plaquette lattices. The plaquette lattice has a square plaquette in each unit cell 
with four lattice-points at the vertices, as shown in Fig. 1(b). It can be shown that both 
QDSLs are well represented by the Hubbard model. An interesting difference between the 
two QDSLs is that the Fermi surface of the plaquette QDSL has disconnected pieces whereas 
the square QDSL is formed of one piece. To find a correlation effect that reflects both the 
Coulomb interaction and the structures of the Fermi surface, we studied the existence of 
superconductivity for both lattices within the framework of the Hubbard model. We found a 
superconducting ground state where the transition temperature T c of the plaquette lattice is 
more than double that of the square lattice and sufficiently high to allow superconductivity 
to be observed experimentally. 

We begin with the design of the square QD lattice. We consider the quantum wire 
network shown in Fig. 1(a). We assume InAs quantum wires buried in In .776Ga .224As 
barrier regions with a band offset of 0.17 eV ||. Each quantum wire is 50 nm wide and 
50 nm high, and the lateral distance between adjacent wires is 61.1 nm. These parameters 
have been carefully determined. If the lateral distance is too small, the wave function of an 
electron is extended and QDs are not appropriately formed. By contrast, if the distance is 
too large, the wave function is too localized at the intersection and the hopping parameter 
between adjacent QDs becomes very small. The electronic band structures are obtained 
using first-principles calculations based on the LDA with the Perdew-Zunger exchange- 
correlation potential [10]. We assume the effective mass of electrons to be m* = 0.02m and 



the dielectric constant e = 12.4eo for the InAs wires. We employ conventional plane wave 
expansion for the LDA calculation, and include a sufficiently large number of plane waves 
to ensure numerical convergence. 

The LDA band diagram for this wire network is shown in Fig. 2. It can be perfectly 
fitted to the tight-binding calculation for the square lattice parametrized in Fig. 1(b): the 
fitting parameters are hopping integrals between adjacent sites {t a = U = 0.179 meV) and 
much smaller ones between non-adjacent sites (t c = = tf = —0.017 meV and t e = —0.021 
meV). This fitting result indicates that the tight-binding approximation is very accurate and 



the wave function is appropriately localized on a QD ||11|| . The localization originates from 



the coherence between the crossed plane waves at the intersection, resulting in a localized 
state with reduced energy. From another point of view, because of the uncertainty condition 
between position and momentum in quantum mechanics, the momentum of electrons, or the 
kinetic energy, becomes lower at a wider area in the intersection. Note that the results for 
fitted hopping parameters between the nearest neighbor QDs are sufficiently large and the 
wave function is not too localized on a QD. We should also note that the fitted hopping 
parameters between the second nearest dots (t c , td, and £/) and those between the third 
nearest dots (t e ) are small. This is due to the large band-offset of Ino.776Gao.224As which 
suppresses tunneling through this region. Therefore, the hopping parameters are determined 
by the distance not in a straight line but along the wires. The localization of the wave 
function also suggests that the intra-dot Coulomb interaction between electrons is large. 
The intra-dot interaction, simply estimated by subtracting the Ewald sum for inter-dot 
interactions from the total Hartree energy in the LDA, is U = 1.8 meV, which is similar to 
that estimated in a QD formed in a spherically symmetric harmonic potential with a 50 nm 
oscillator length. The value of U/t a ~ 10 shows that the electrons are strongly correlated 
and the many-body effects of the electrons should be pronounced. 

Next, we consider the plaquette QDSL [Fig. 1(a)], which we obtained by alternating 
the distance between the wires in the wire network for the square QDSL. In our design, 
the smaller (larger) distance between adjacent wires is 38.8 (83.4) nm. In the tight-binding 
model parametrized in Fig. 1(b), the hopping parameter between adjacent sites within a 
plaquette and the one across two plaquettes are different (t a ^ tb) unlike in the square lattice 
{t a = tb)- By performing the same analysis undertaken for the square lattice, we confirmed 
that the tight-binding model provides a good fit with the LDA band diagram with the fitting 
parameters t a = 0.242, t b = 0.151, t c = 0.003, t d = -0.025, t e = -0.024, and t f = -0.014 in 
milli-electron volt units. The estimated intra-dot interaction is 1.7 meV. 

A characteristic difference between the square QDSL and the plaquette QDSL lies in the 
structure of the Fermi surface shown in Fig. 3, which is calculated in the LDA with the elec- 



tron filling n = 0.8 (=number of electrons/number of QDs). The Fermi surface of the square 
QDSL is connected, while that of the plaquette QDSL has disconnected pieces originating 
from the band folding. To demonstrate that the difference in the structure of the Fermi 
surface significantly influences the characteristics of the correlated electrons in QDSLs, we 
studied superconducting correlations in both QDSLs by assuming the Hubbard model with 
the estimated hopping and intra-site interaction parameters described above. It has been 
shown that the repulsive Hubbard model exhibits a superconducting transition. In addition 
to early calculations |12| , a recent quantum Monte Carlo calculation indeed shows the 



enhancement of the pairing correlation with <i-wave symmetry in the two-dimensional re- 
pulsive Hubbard model. Here, we employ the fluctuation-exchange (FLEX) approximation 
|14Hl9| with the Eliashberg equation [|20| in order to evaluate the superconducting transi- 
tion temperature T c in the Hubbard model for the plaquette and square lattices. The FLEX 
approximation is a kind of self-consistent random phase approximation and is known to be 
appropriate for treating the strong antiferromagnetic spin fluctuation that originates from 
the Coulomb interaction between electrons with opposite spin. We note that the FLEX ap- 
proximation gives a good estimate of T c for cuprates [P^T8| . For both lattices, we assume 



the filling of electrons n = 0.8, where the system is far from the antiferromagnetic instabil- 
ity region near half-filling (n = 1) but the spin fluctuation is still large enough to obtain 
a high T c . Our analysis shows that T c for the plaquette QDSL is 90 mK, which is more 
than double that of T c =40 mK for the square QDSL. Both T c values are high enough for 
us to observe superconductivity experimentally. However, low-temperature measurements 
are usually limited to a few tens of milli-Kelvins, and a higher T c would be preferable as 
it would allow us to observe the superconducting transition more clearly. Figure 4 shows 
the improvement in T c when the hopping parameters t a and % are changed for n = 0.85 
and U = 7. Other small hopping parameters are disregarded {t c = t d = t e = tj = 0) 
because they have little influence on the result. Energy is scaled in units of average hopping 
energy t — (t a + t b )/2 = 1. In these units, the single-particle bandwidth is simply given by 
4(t a + h) = 8. We find that T c is enhanced as the difference between t a and U is increased, 



and T c for t a = 1.5 and U = 0.5 is more than double that for the square lattice (t a = U = 1). 

The mechanism of this large T c enhancement is similar to that recently proposed by 
Kuroki and Arita |2T] . The strength of the superconducting pairing interaction is charac- 
terized by the coupling constant in the form 

Z)fc,fc' e Fermi surfaces H>air(& ~ k)(j){k)(j){k') 
►efT — ~ [ -L ) 

6(k) 2 

^fcGFermi surfaces t\ ) 

where V va \ r is the effective electron-electron interaction induced by the pair-scattering pro- 
cess from [k, —k] to [£?, —k'j on the Fermi surface and <p(k) is the superconducting gap 
function. Positive values of V e R correspond to an attractive interaction and, roughly speak- 
ing, a larger V e s gives a higher T c . Reflecting the Fermi surface nesting in the Hubbard 
model, antiferromagnetic spin fluctuations near half filling enhance the spin susceptibility 
x(q) [roughly proportional to Vp air (g*)] around a certain wave vector q = Q. When pair 
scattering on the Fermi surface is accompanied by a momentum transfer k — k! around the 
peak Q and a sign change of the gap function cj)(k)(j)(k') < 0, this scattering process gives 
significant positive contributions to V e g. In the square lattice, the Fermi surface is perfectly 
nested at half-filling and x(o) exhibits the antiferromagnetic peak at Q = Since the 

gap function <p{k) = cos k x — cos k y with d x 2_ y 2 symmetry has nodes on the Fermi surface, 
some pair scatterings transferring k — k' = Q w (7r, it) between two points [indicated by solid 
arrows in Fig. 3(a)], where the gap function has the opposite sign, make significant positive 
contributions to V e s. However, some other scatterings k — k' — Q + AQ with the same sign 
of 4>{k) (indicated by dotted arrows) reduce V e Q. In the plaquette lattice, the Fermi surface 
is folded in half and an antiferromagnetic peak emerges at Q ~ with a finite spread AQ. 
The Fermi surface is separated by lines \k x \ = \k y \, and is separated into pieces, each of 
which has a gap function with the same sign, as shown in Fig. 3(b). Dominant scattering 
processes around Q w always make positive contributions to Kff and significantly enhance 
T c . For a dimerized lattice, which is somewhat similar to a ladder system |[22|| , Kuroki and 
Arita pi[ have shown that a very high T c can be reached based on this mechanism. However, 



to raise T c in a dimerized lattice, it is necessary to tune three different hopping-parameters, 



whereas, in plaquette lattice-patterns, only two hopping parameters determine T c , and these 
can be easily modulated by selecting the wire geometry. 

In summary, we proposed a QDSL design based on square and plaquette type structures 
in quantum wire networks. Both QDSLs are well represented by the Hubbard model. We 
studied superconducting transitions for both QDSLs using the Hubbard model and found 
the T c for the plaquette QDSL to be more than double that for the square QDSL. The 
estimated T c value in the plaquette QDSL is sufficiently high to allow us to observe super- 
conductivity experimentally. Finally, we note that the transition temperature of cuprates 
might be increased close to room temperature if the material is patterned in a plaquette 
type lattice. 
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Aoki for stimulating discussions. They also thank A. Richter for his critical reading of the 
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FIGURES 

FIG. 1. (a) Diagram of a nonuniform quantum wire network with distances a and b between 
adjacent wires, (b) Corresponding tight-binding model with hopping parameters t. For a square 
QDSL, t a = tb and t c = t d = tf. 

FIG. 2. Result of fitting the band structures of the square QDSL. The LDA band diagram 
of a quantum wire network (solid curve) and the best fit in the tight-binding calculation (dashed 
line) with parameters t a = tb = 0.179 meV, t c = td = tf = —0.017 meV, and t e = —0.021 meV are 
shown. 

FIG. 3. Energy contours for (a) the square QDSL and (b) the plaquette QDSL calculated in 
LDA. The energies are plotted in milli-electron volts. The thick solid (dashed) curves represent 
the Fermi surfaces where the superconducting gap function is positive (negative). The upward 
(downward) pointing arrow indicates an electron with up (down) spin. The solid (dashed) arrow 
schematically represents dominant pair-scattering processes that give a positive (negative) contri- 
bution to the pairing interaction by transferring a momentum at which the spin susceptibility has 
a peak with a finite spread. 

FIG. 4. T c plotted as a function of t a in units of the average hopping parameter 
t = (t a + *b)/2 = 1. The single-particle bandwidth is given by 4(t a + tb) = 8. Parameters 
t a = h = 1 (f° r the square lattice) gives T c = 0.025. 
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